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Abstract 

Solar flares are accessible to a broad variety of observational methods to see and inves- 
tigate the magnetic reconnection phenomenon in high-temperature strongly-magnetized 
plasma of the solar corona. An analysis of the topological peculiarities of magnetic field 
in active regions shows that the topological trigger effect is necessary to allow for in order 
to construct models for large eruptive flares. 

The topological trigger is not a resistive instability which leads to a change of the 
topology of the field configuration from pre- to post reconnection state. On the contrary, 
the topological trigger is a quick change of the global topology, which dictates the fast 
reconnection of collisional or collisionless nature. 

The current state of the art and development potential of the theory of collision- 
less reconnection in the strong magnetic fields related to large flares are briefly reviewed. 
Particle acceleration is considered in collapsing magnetic traps created by reconnection. 
In order to explain the formation of coronal X-ray sources, the Fermi acceleration and 
betatron mechanism are simultaneously taken into account analytically in a collisionless 
approximation. Finally, the emphasis is on urgent unsolved problems of solar flare physics. 

1 The role of magnetic fields in flares 
1.1 Basic questions 

Frequent observations of solar flares became available in the 1920s. Early studies showed 
that flares were associated with magnetic fields. Estimates of the energy required to power 
large flares led to the conclusion that flares must be electromagnetic in origin. Step by step 
it became more clear that a flare is the result of reconnection of magnetic fields. However 
there were some objections to the hypothesis that the energy of a flare could be stored in 
the form of a magnetic field of reconnecting current layers (RCLs). 

(1) It was claimed that measurements of photospheric magnetic fields do not demon- 
strate an unambiguous relation between flares and the changes of the fields. More exactly, 
the changes in question are those that occur immediately before a flare to create it. These 
changes were supposed to be the cause of a flare. 
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(2) Next objection was related to the time of dissipation of the magnetic field in a 
volume that would contain the energy necessary for a flare. If this time is estimated in 
a usual way as the diffusion time in a solar plasma of a finite conductivity, then it is too 
long compared with the observed duration of a flare. 

(3) The third objection was the most crucial one: the observers had never seen real 
RCLs in solar flares. 

Starting from Severny [1], solar observers have been studying flare-related changes in 
photospheric magnetic fields to provide information how an active region (AR) stores and 
releases its energy (see Lin et al. [2], Wang [3]). However the role of photospheric fields is 
still an area of ongoing research (Section ll.3p . What are the answers of the reconnection 
theory to the objections listed above? 

1.2 Concept of magnetic reconnection in solar flares 

According to contemporary views, the principal flare process is contingent on the accu- 
mulation of the free magnetic energy in the corona. By 'free' we mean the surplus energy 
above that of a potential magnetic field 



Here ?/>(r) is the potential of the field, the index a = 1,2,3. The potential field has the 
sources (sunspots, background fields) in the photosphere. 

The free magnetic energy is related to the electric currents in the corona. A solar flare 
corresponds to rapid changes of these currents. So we distinguish between two processes: 
a slow accumulation of energy and its fast release, a flare. Let us see them in a classi- 
cal example, the evolution of the quadrupole configuration of sunspots (Figure [I]). The 
sunspots of pairwise opposite polarity are shown: N and S represent a bipolar group of 
sunspots in an AR, n and s model a new emerging flux. All sunspots are placed along the 
axis x in the photospheric plane Ph at the bottom of the chromosphere Ch. 

The field line A\ in Figure QJi is the separatrix line of the initial state (a), this line 
will reconnect first. X is the zeroth point of the field at the initial state, here the RCL is 
created at the state (b). The field line A<i is the separatrix of the final state (c) or the last 
reconnected field line. Therefore 5 A = A2 — A\ is the reconnected magnetic flux. 

Three solid arrows under the photosphere Ph in Figure [lb show a slow emergence of 
a new magnetic flux (the sunspots n and s). The sunspots have been emerged but the 
field lines do not start to reconnect. More exactly, they reconnect too slowly because of 
high conductivity of plasma in the corona. In the first approximation, we neglect this 
reconnection. 

In general, the redistribution of fluxes appears as a result of the slow changes of the 
field sources in the photosphere. These changes can be either the emergence of a new flux 
from below the photosphere (Figure [1]) or other flows of photospheric plasma, in particular 
the shear flows, inhomogeneous horizontal flows along the neutral line of the photospheric 
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Figure 1: The classical 2D cartoon of magnetic reconnection in a solar flare. Three states 
of the field: (a) the initial one, (b) the pre-reconnection state, (c) the final state after 
reconnection. 

magnetic field. For this reason, an actual (slow and fast) reconnection in the corona is 
always a three-dimensional (3D) phenomenon (Section [2]). 

1.3 Some new results of the magnetic field observations 

Let us come back to objection (1) in Section 11.11 to the reconnection theory of flares. 
According to the theory, the free energy is related to the current J inside the RCL. A 
flare corresponds to rapid changes of this current. However the magnetic flux through the 
photosphere Ph (Figure [1]) changes only little over the whole area of a flare during this 
process, except in some particular places, for example, between close sunspots N and s. 
Sunspots in the photosphere are weakly affected by a flare because the plasma in the 
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photosphere is almost 10 9 times denser than the plasma in the corona. It is difficult 
for disturbances in the tenuous corona to affect the extremely massive plasma in the 
photosphere. Only small perturbations penetrate into the photosphere. 

The same is true in particular for the vertical component of the field. The photospheric 
magnetic field changes a little during a flare over its whole area. As a consequence, 
after a flare the large-scale structure in the corona can remain free of noticeable changes, 
because it is determined mainly by the potential part of the field above the photosphere. 
More exactly, even being disrupted, the structure will come to the potential configuration 
corresponding to the post-flare position of the photospheric sources. 

On the other hand, in the "Bastille-day" flare on 2000 July 14 and in some other large 
flares, it was possible to detect the real changes in a sunspot structure just after a flare. 
The outer penumber fields became more vertical due to reconnection in the corona during 
a flare (Liu et al. [4], Wang et al. [5]). One can easily imagine such changes by considering 
Figure [1] between sunspots N and s. 

Sudol and Harvey [6] used the magnetograms to characterize the changes in the ver- 
tical component of photospheric field during 15 large flares. An abrupt, significant, and 
persistent change occurred in at least one location within the flaring AR during each event 
after its start. Among possible interpretations, Sudoh and Harvey favor one in which the 
field changes result from the penumber field relaxing upward by reconnecting magnetic 
field above the photosphere. This interpretation is similar to than one given by Liu et 
al. [4] and Wang et al. [5]. 

As for objection (2) to the hypothesis of accumulation of energy in the form of magnetic 
field of current layers, the rapid dissipation of the field necessary for a flare is explained 
by the theory of super-hot turbulent-current layers (Section [4]) . In what follows, we advo- 
cate that the potential part of magnetic field in the solar corona determines a large-scale 
structure and properties of active regions, while the reconnecting current layers determine 
energetics and dynamics of flares. 

2 Three-dimensional reconnection in flares 

2.1 The first topological model of an active region 
2.1.1 The potential field approximation 

Gorbachev and Somov [7,8] have developed a 3D model for a potential field in the AR 2776 
with an extended flare of 1980 November 5. Before discussing the AR and the flare (see 
Section I2T21) . let us consider, at first, the general properties of this class of models called 
topological. 

In the simplest model (Gorbachev and Somov [9]), four field sources - the magnetic 
"charges" e N and e s , e n and e s , located in the plane Q under the photosphere Ph (Figure[2]) 
- are used to reproduce the main features of the observed field in the photosphere related 
to the four most important sunspots: N, S, n and s. As a consequence, the quadrupole 
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Figure 2: A model for the magnetic field of four sunspots of pairwise opposite polarity. 
The sunspots iV and S in the photospheric plane Ph. The separatrices Si and 52 cross 
at the separator X1XX2 above the plane Q of the effective magnetic charges e N , e s , etc. 

model reproduces only the large-scale features of the actual field in the corona related to 
these sunspots. As a minimum, the four sources are necessarily to describe two interacting 
magnetic fluxes having the two sources per each. The larger number of sources are not 
necessarily much better. 

The main features are two magnetic surfaces called the separatrices: S\ and S2 (Fig- 
ure^. They divide the whole space above the plane Q into four regions and, correspond- 
ingly, the whole field into four magnetic fluxes having different linkages. The field lines 
are grouped into four regions according to their termini. The separatrices are formed from 
lines beginning or ending at magnetic zeroth points X\ and X^. For example, the field 
lines originating at the point X% form a separatrix surface S\. 

The topologically singular field line X1XX2, lying at the intersection of the separat- 
rices, belongs to all four fluxes (two reconnecting and two reconnected fluxes) that interact 
at this line, the 3D magnetic separator. So the separator separates the interacting fluxes 
by the separatrices. Detection of a separator, as a field line that connects two zeroth 
points in the Earth magnetotail, by the four Cluster spacecraft (Xiao et al. [10]) provides 
an important step towards establishing an observational framework of 3D reconnection. 

2.1.2 Classification of zeroth points 

Let us clarify properties of the zeroth points of the magnetic field discussed above. In 
the vicinity of a zeroth point Xi located at r = rj in the plane Q, the vector of field is 
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represented in the form 

B a (r) = J^- =M aP (n)^. (2) 

Here tp( r ) is the potential of the field, the vector £ = r — T{, the index i = 1, 2 ; the Greek 
indices a, /3 = 1,2,3. The symmetric matrix 



dr Q drp 



Ai 
A 2 
A 3 



(3) 



in the frame of coordinates r' a related to the eigenvectors e^, X a are the eigenvalues of the 
matrix. Because the field is potential, all the \ a are real numbers. Positive eigenvalues A a 
correspond to the field lines emerging from a zeroth point and negative X a correspond to 
the lines arriving at a point rj. 
If the determinant 

detM a/3 (r i ) = A 1 A 2 A 3 /0, (4) 

then the zeroth point is called non- degenerate. Such a point is isolated, i.e., the magnetic 
field in its vicinity does not vanish. Since div B = 0, we obtain a condition 

Ai + A 2 + A 3 = 0. (5) 

Therefore there exist only three following classes of zeroth points. 

Type A. All eigenvalues A Q ^ 0. This corresponds to two zeroth points shown in 
Figure El At the first one X \ 

Ai<0, A 2 >0, A 3 >0. (6) 

Using the analogy with fluid flow, we say that field lines "flow in" along the x axis and 
"flow out" along the separatrix plane (y, z) as shown in the left bottom insert in Figure [2j 
We shall call this subclass of non-degenerate zeroth point as the Type A— since A i < . 
On the contrary, at the point AT 2 

Ai>0, A 2 <0, A 3 <0. (7) 

This will be the Type A+ since A i > . 
Type B. This is a degenerate case 

Ai = -A 2 /0, A 3 = 0. (8) 

It can occur, for example, if a zeroth point belongs to a neutral line, i.e. a line consisting 
of zeroth points of the hyperbolic type. 

Type C. This is another degenerate case 

Ai = A 2 = A 3 = 0. (9) 

It means that three neutral lines intersect at the point rj . 
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Table 1: Topological indices of singular and zeroth points 



Type of a point 


Itop 


Maximum of potential tp 


+1 


Minimum of potential ip 


-1 


Zeroth point of Type A+ 


+1 


Zeroth point of Type A— 


-1 



2.1.3 The number of zeroth points 

In order to use some general theorems of differential geometry, we have to smooth the point 
charges of magnetic field. The smoothed potential ip of a positive source has a maximum 
at the point where this charge is located. Hence at this point all eigenvalues A Q > 0. At 
a singular point with a negative charge, all X a < 0, and the potential has a minimum. 

In order to establish a relation between the number of non-degenerate points, let us 
introduce the topological index 

ho P (?i) = sign[detM a/3 (ri)] = sign(A x A 2 A 3 ) . (10) 

Here the function sign (x) = +1 if x > and sign (x) = —1 if x < 0. Thus the possible 
types of points have the indices presented in Table HJ 

According to Dubrovin et al. [11] (see Part II, Ch. 3, § 14), with some difference in 
notations that are standard in physics, we formulate the following statement. 

Theorem 1: For a 3D potential magnetic field, in a general 3D case of magnetic source 
location, 

V It (r ) - J- [B x ™1 - J® (11) 

^ Itop[ri} ~ Att J ^ [ d9 X dip\ ~ J ' (Uj 

where the integral is taken over the sphere of infinite radius. 
If the total magnetic charge 

e to4 = 5> 4 >0, (12) 

i 

then the integral J® = +1. Let N max (N m i n ) is the number of maxima (minima) of the 
potential ip of magnetic field B, (Na~) is the number of zeroth points of Type A+ 

(A—), then according to (jlip we have 

N max - N min + N A+ - N A _ = +1 . (13) 
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If etot < then = — 1, and we write —1 on the right-hand side of Equation (fT3|) , 
If the total charge 

e t ot = 0, (14) 
then J ® = 0, and the following relation is valid: 

N max - N min + N A+ -N A „ = 0. (15) 

For the magnetic field shown in Figure [21 N max = N m i n = 2. Hence formula ([15ft gives 
us an equation 

N A+ = iV A _ (16) 

but not the total number of zeroth points, that we need. 

However, if all the charges are located in the plane Q as illustrated by Figure [21 then we 
can obtain an additional information about the number of zeroth points, another equation. 
Let us determine the topological index for the 2D vector field in the same plane 

B®(x,y) = {B x ,B y } (17) 
following general definition from Dubrovin et al. [11] 

l/ 5(r i ) = sign(A 1 A 2 ) . (18) 
Hence the possible types of non-degenerate points have the indices presented in Table [21 

Table 2: Topological indices of singular and zeroth points of the 2D field in the plane Q 



Type of a point 


r (2) 

1 top 


Maximum of potential ip 


+1 


Minimum of potential tp 


+1 


Zeroth point (a saddle) 


-1 



In this case of a plane arrangement of magnetic charges, we have to rewrite Theorem 1 
as follows. 

Theorem 2: For a potential magnetic field in the plane Q, in which the charges are 
located, 



(2) 
top 1 



1 

2^ 



B ( 2 ) {x, y) 2 ( B x dB y - B y dB x ) = J 



(2) 



(19) 



Here the integral is taken over the circle of infinite radius in the positive direction of 
circulation. 

For a plane field (|T7|) of the dipole type (i.e., etot = and an effective dipole mo- 
ment m/0) the integral = 2. Thus we have 

Nit + Ni 2 l-N[ 2) = 2, (20) 

where Nrkax {N m - n ) is the number of maxima (minima) of the potential ip = ip{x,y,0) in 

(2) 

the plane z = 0, N A is the number of zeroth points, the saddles. Because of the symmetry 
relative to the plane Q, everywhere in the plane (outside of the charges) B z = 0. Therefore 
the saddles of the plane field are zeroth points of a 3D field. 

(2) (2) 

For the field shown in Figure [21 N m dx = 2 and N^' n = 2. Hence, according to 

formula (|20p . the number of zeroth points N\ =2. By using Equation (|16p . we conclude 
that, in the simple case under consideration, 

N A+ =N A - = 1. Q.e.d. (21) 



2.1.4 The electric currents needed 

The potential field model does not include any currents and so cannot model an energy 
stored in the fields and released in flares. Therefore here we introduce some currents and 
energetics to a flare model. 

The topological model under consideration does not mean, of course, that we assume 
the existence of magnetic charges under the photosphere as well as the real zeroth points X\ 
and X2 in the plane Q which does not exist either. We assume that above the photosphere 
the large-scale field can be described in terms of such a simple model. If the magnetic 
sources move and change, the field also changes. It is across the separator that the 
magnetic fluxes are reconnected so that the field could remain potential, if there were no 
plasma. 

In the presence of a plasma of low resistivity, the separator plays the same role as the 
hyperbolic zeroth line of magnetic field, familiar from 2D MHD problems (Syrovatskii [12], 
Brushlinskii et al. [13], Biskamp [14,15]; see review of a current state of numerical sim- 
ulations and laboratory experiments in Yamada et al. [16]. In particular, as soon as the 
separator appears, the electric field Eo induced by the varying magnetic field produces 
an electric current J along the separator. The current interacts with the potential field 
in such a way that the current assumes the shape of a thin wide current layer (RCL in 
Figure EJ. 

In a low-resistivity plasma the current layer hinders the redistribution of magnetic 
fluxes, their reconnection. This results in an energy being stored in the form of magnetic 
energy of a current layer - the free energy. Therefore the slowly-reconnecting current 
layer appears at the separator (Syrovatskii [17], Longcope and Cowley [18]) in a pre-flare 
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Figure 3: A schematic view of a RCL with a total current J at the separator. 



stage. If for some reason (see Somov [19], Cassak et al. [20], Uzdensky [21]) reconnection 
becomes fast, the free magnetic energy is rapidly converted into energy of particles. This 
is a flare. The rapidly-reconnecting current layer, being in a high-temperature turbulent- 
current state (Section S]), provides the energy fluxes along the reconnected field lines. 

Reconnection for 3D models without zeroth points in and above the photosphere, as we 
assumed in Figure is qualitatively the same as the 2D reconnection in the models with 
a non-zero longitudinal magnetic field B n . The inductive electric field along the separator 
is the driving force of reconnection. The field lines going to reconnect are pushed into the 
RCL and the lines coming out of the RCL are pulled by the longitudinal electric field. 
Thus reconnection in 3D should be taken to mean effects associated with the RCL, as in 
2D (Greene [22], Lau and Finn [23]). 

An actual 3D reconnection at a separator proceeds in the presence of an increasing (or 
decreasing) longitudinal field By. What factors do determine the increase (or decrease) 
of this field? - The first of them is the global field configuration, i.e. the relative position 
of the magnetic field sources in an AR. It determines the position of a separator and the 
value of the longitudinal field at the separator and in its vicinity. This field is not uniform, 
of course. 

The second factor is evolution of the global configuration, more exactly, the electric 
field Eo related to the evolution and responsible for driven reconnection at a separator. 
The direction of reconnection - with an increase (or decrease) of the longitudinal magnetic 
field - depends on the sign of the electric field projection on the separator, i.e. on the sign 
of the scalar product (Eo • Bn ). In general, this sign can be plus or minus with equal 
probabilities, if there are no preferential configurations of the global field or no preferential 
directions of the AR evolution. 
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2.2 The solar flare on 5 November 1980 
2.2.1 Observed and model magnet ograms 

The first well-studied example was the extended 1B/M4 flare on 5 November 1980 (e.g. 
Rust and Somov [24] ) observed by the satellite SMM (Figure 0|) . Three bright hard X- 
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Figure 4: The flare of 1980 November 5 as observed by SMM . 

ray (HXR) kernels A, B and C are well distinguished in the right panel. The soft X-ray 
(SXR) elements of the flare, as seen in the left panel, consist presumably of two overlapping 
coronal loops, AB and BC. Their footpoints in the chromosphere have been identified as 
the sources of HXR and Ha emission. 

Figure [5] shows the line-of-sight magnetogram of the AR 2776 where the flare occurred. 
Numbers shows a magnitude of the field measured in 10 2 G. Two narrow flare ribbons 
are shown on either side of the photospheric neutral line (PNL) . Let us identify the four 
largest regions in which the field of a single polarity is concentrated: two of northern 
polarity and two of souther polarity. Since the AR is comparatively close to the center of 
the solar disk, the magnetogram represents the vertical component B z at the photospheric 
level fairly well. 

Let us make a model of the magnetogram shown in Figure O With this aim, we calcu- 
late the vertical component B z of magnetic field in the photospheric plane Ph (plane z = 
in this Section) in the potential approximation by using the simple formula 

B (*, y,z)=J2 | e% ,2 x T~~^l • (22) 
1=1 1 1 1 1 1 1 

Here are the effective charges of field, and rj are their radius vectors: 

{-0.8,0 ,-0.1}; 
{0.1,-0.4,-0.1}; 
{-0.25,0.15,-0.1}; 
{0.5,0.16,-0.1} . 



ei 


= -2.5, 


r-i 


C J 2 


= -1.0, 






= +3.0, 




e 4 


= +4.5, 


ri 
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Figure 5: The line-of-sight component of magnetic field in AR 2776. Two flare ribbons 
observed in Ha line are shown as shadow areas; FR\ is the longer one. 

Note that the total magnetic charge 



Thus, in order to determine a relation between the number of zeroth points and the 
numbers of field sources, we have to use formula (fT3|) . It gives us an equation 



We shall come back to this equation, considering a topological portrait of the AR. 

The parameters e$ and r-j have been selected in order to reproduce in an optimal way 
the fluxes of individual polarities in the photospheric plane (Figure [6^) as well the shape 
of the photospheric neutral line (PNL in Figure ©)) . The calculated simplified neutral 
line (SNL) smoothes the curve PNL in small scales but conserves its S"-type shape in 
a large scale, in the linear scale of the AR. In Figured! the length unit equals 150" ~ 
1.1 x 10 10 cm; the numbers near the curves show the values of the vertical component of 
the field measured in 10 2 G. 

Because the simplest model uses a minimal number of magnetic sources - four, which 
is necessary to describe the minimal number of interacting fluxes - two, we call it the 
quadrupole-type model (Figure [6]). This is not an exact definition (because in the case 
under consideration e N ^ — e s and e n / — e s ) but it is convenient for people who know 
the exact- quadrupole model (Sweet [25]). In fact, the difference, the presence of another 



e to t = 4 > 0. 



(23) 



N A+ = iV A _ + 1 . 



(24) 
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Figure 6: The model magnetogram. (a) Positions and magnitudes of the sources fit 
the main features of the observed magnetogram. (b) The calculated simplified neutral 
line (SNL) approximates the observed photospheric neutral line (PNL). 

separator in the model under consideration, is not small and can be significant for actual 
ARs. The second separator may be important to give accelerated particles a way to escape 
out of an AR in the interplanetary space. 

We calculate the field lines integrating the ordinary differential equations 



Here (dx, dy, dz) is an arch element directed along a field line. The vector B is determined 
at each point by formula (|22|) . 

2.2.2 Topological portrait of the active region 

The topological properties of magnetic field are determined by the number and locations 
of zeroth points. In the case under consideration, there are three points in the plane Q 
(here z = —0.1) shown in Figure [3 The coordinates of these points are: 



Figure [7] shows the topologically important field lines in the plane (x',y') which is 
the plane Q of the sources ei, 62, e3, and e^. Since the total charge etot > 0, instead of 



dx dy dz 

B X By B z 



(25) 



X x = {-0.133,-0.739,-0.1} 
X 2 = {0.076, 0.179,-0.1}, 
A 3 = {-2.26,-0.439,-0.1}. 
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Figure 7: The topological portrait of the AR 2776 where the extended 1B/M4 solar flare 
of 1980 November 5 occurred. 

Equation (|20p for the number of zeroth points of 2D field (117H . we have an equation 

NH + N^l-Ni 2) = 1, (26) 

because according to formula (]19p we find the value J (2) = 1. It follows from (J26D that 
the total number of zeroth points 

ivj 2) = 3 . (27) 

On the other hand, 

ivf = N A+ + N A - . (28) 

With account of ([23]) taken, we conclude that Na+ = 2 and Na~ = 1. 

The field lines shown in Figure [7] play the role of separatrices and show the presence 
of two separators. Two points X\ (Type Af) and X2 (Type A_) are located in the 
vicinity of the magnetic sources and are connected by the first separator shown by its 
projection, the thin dashed line L\. Near this main separator, the field and its gradient 
are strong and determine the flare activity of the AR. 

Another separator starts from the point A3 (the Type A + ) far away from the sources 
and goes much higher above the AR to infinity. The second separator can be responsible 
for flares in weaker magnetic fields and smaller gradients high in the corona. The second 
separator is also a good place for the particles accelerated along the main separator to 
escape from the AR in the interplanetary space. 
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Figure 8: Potential field lines cross a region of primary energy release £ situated at the 
apex of the main separator (boldface dashed curve). Two flare ribbons are formed where 
these lines cross the photospheric plane Ph. 

We suppose that a part of the flare energy is released in some region £ near the apex 
of the main separator. The energy fluxes F £ propagate along the field lines connecting the 
energy source with the photosphere. Projections of the energy source £ on the photospheric 
plane Ph along the field lines are shown as two "flare ribbons" FR\ and FR2 in Figure [8l 
Therefore the model allows us to identify the flare brightenings, in the Ha line as well as 
in EUV and HXRs, with the ribbons located at the intersection of the separatrices with 
the chromosphere which is placed slightly above the photospheric plane (Ph). 

The characteristic saddle structure of the orthogonal (perpendicular to a separator) 
field Bj_ in the vicinity of the separator leads to a spatial redistribution of the energy 
flux F £ of heat and accelerated particles. This flux is split apart in such a way that it 
creates the long-narrow ribbons. Thus, at first, the model had reproduced the observed 
features of the real solar flare. In particular, the model predicts the simultaneous flaring of 
the two chromospheric ribbons. Moreover it predicts that a concentration of the field lines 
that bring energy into the ribbons is higher at the edges of the ribbons, i.e. at relatively 
compact regions indicated as A, B, and C. Here the Ha brightenings must be especially 
bright. This is consistent with observations of Ha "kernels" in this flare. 

2.3 Features of the flare topological model 

Since in the kernels the energy fluxes are concentrated, the impulsive heating of the chro- 
mosphere creates a fast expansion of high-temperature plasma upwards into the corona 
(e.g. Somov [26]). This effect is known as the chromospheric "evaporation" observed in 
the EUV and SXR emission. Evaporation lights up the coronal loops in flares. 

The topological model also shows that the flare ribbons as well as their edges with 
Ha kernels are connected to a common source of energy at the separator (see £ in Figure [9]). 
Through this region all four kernels are magnetically connected to one another. Therefore 
the SXR loops look like they are crossing or touching each other. 
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Figure 9: Field lines that connect the Ha kernels A, B, C, and D. Chrome-spheric 
evaporation creates a picture of the crossing SXR loops. 

So the reconnecting magnetic fluxes are distributed in the corona in such a way that 
the two SXR loops may look like that they interact with each other. That is why 
such structures are usually considered as evidence in favor of a model of two interacting 
loops (Sakai and de Jager [27]). The difference, however, exists in the primary source 
of energy. High concentrations of electric currents and twisted magnetic fields inside the 
interacting loops are created by some under-photospheric dynamo mechanism. If these 
currents are mostly parallel they attract each other giving an energy to a flare (Gold and 
Hoyle [28]). On the contrary, according to the topological model, the flare energy comes 
from an interaction of magnetic fluxes that can be mostly potential. 

The S-shaped structures observed in SXRs are usually interpreted in favour of 
non-potential fields. In general, the shapes of loops are signature of the helicity of their 
magnetic fields. The S-shaped loops match flux tubes of positive helicity, and inverse 
S-shaped loops match flux tubes of negative helicity (Pevtsov et al. [29]). However, in 
Figure [H the S-shaped structure CEB connecting the bright points C and B results from 
the computation of the potential field. 

In the AR 2776, Den and Somov [30] found a considerable shear of a potential field 
above the PNL near the brightest loop AB. Many authors concluded that an initial energy 
of flares is stored in magnetic fields with large shear. However, such flares presumably were 
not the case of potential field having a minimum energy. This means that the presence of 
observed shear is not a sufficient condition for generation of a flare. 

The topological model by Gorbachev and Somov [7] postulated a global topology for 
an AR consisting of four main fluxes. Reconnection between, for example, the upper and 
lower fluxes transfers a part of the magnetic flux to the two side systems. Antiochos [31] 
addresses the following question: What is the minimum complexity needed in the magnetic 
field of an AR so that a similar process can occur in a fully 3D geometry? He starts with a 
highly sheared field near the PNL held down by an overlying unsheared field. Antiochos 
concludes that a real AR can have much more complexity than very simple configurations. 
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We expect that the large-scale topology of four-flux systems meeting along a separator is 
the basic topology underlying eruptive activity of the Sun. 

On the other hand, configurations with more separators (as well as with multiply 
connected magnetic sources and zeroth-point pairs; see Parnell [32]) have more opportunity 
to reconnect and would thus more likely to produce flares. Such complicated configuration 
would presumably produce many small flares to release a large excess of energy in an AR 
rather than one large flare. 

3 Topological trigger for solar flares 
3.1 What is that? 

The effect of topological trigger was suggested by Gorbachev et al. [33] in spite of a 
large language barrier between its mathematical background (Dubrovin et al. [11]) and 
the physical terms used by models of ARs. Many people simply understood the simple 
topological model but, unfortunately, not the topological trigger. 

Fortunately, Barnes [34] investigated a relationship between solar eruptive events and 
the existence of the coronal zeroth points of field, using a collection of over 1800 vector 
magnetograms. Each of them was subjected to the charge topology analysis, including 
determining the presence of coronal zeroth points. It appears that the majority of events 
originate in ARs above which no zeroth point (a magnetic null) was found. However a 
much larger fraction of ARs, for which a coronal zeroth point was found, were the source 
of an eruption than ARs for which no zeroth point was found. Clearly the presence of a 
coronal zeroth point is an indication that an AR is more likely to produce an eruption, 
as 35 % of the ARs for which such a point was found produced eruptions, compared with 
only 13 % of ARs for which no zeroth point was found. We consider this fact as indication 
that the topological trigger can play a significant role in the origin of eruptive flares. This 
is also consistent with the study of Ugarte-Urra et al. [35] . 

The possibility of a topological approach to the question of the trigger for flares was not 
often discussed in the literature. Syrovatskii and Somov [36] considered a slow evolution 
of coronal fields and showed that, during such an evolution, some critical state can be 
reached, and fast dynamic phase of evolution begins and is accompanied by a rapid change 
of magnetic topology. For example, it is possible a rapid 'break-out' of a 'new' magnetic flux 
through the 'old' coronal field of an AR (Figure [TO)) . When an effective magnetic moment 
of the internal growing group of sunspots becomes nearly equal an effective moment of the 
AR (Syrovatskii [37]), the closed configuration quickly turns into the open one. 

In this paper we discuss another possibility. Near a separator the longitudinal com- 
ponent Bu dominates because the orthogonal field vanishes at the separator. Recon- 
nection in the RCL at the separator just conserves the flux of the longitudinal field (see 
Somov [19]). At the separator, the orthogonal components are reconnected. Therefore 
they actively participate in the connectivity change, but the longitudinal field does not. 
Thus it seems that the longitudinal field plays a passive role in the topological aspect of 
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Figure 10: Topological change of a potential field created by four magnetic sources, (a) 
The closed configuration with the zeroth point X. (b) The critical instant when the 
magnetic field opens. 




Figure 11: The zeroth point X c rapidly moves along the separator and switches back the 
longitudinal component By of magnetic field. 



the process but it influences the physical properties of the RCL, in particular the recon- 
nection rate. The longitudinal field decreases compressibility of plasma flowing into the 
RCL. When the longitudinal field vanishes, the plasma becomes "strongly compressible", 
and the RCL collapses, i.e. its width decreases substantially. As a result, the reconnection 
rate increases quickly. However this is not the whole story. 

The important exception constitutes a zeroth point which can appear on the sepa- 
rator above the photosphere. Gorbachev et al. [33] showed that, in this case, even very 
slow changes in the configuration of field sources in the photosphere can lead to a rapid 
migration of such a point along the separator (Figure [TT1) and to a topological trigger 
of a flare. This essentially 3D effect will be considered below in more detail. Note that 
the topological trigger effect is not a resistive instability which leads to a change of the 
topology of the field configuration from pre- to post reconnection state. On the contrary, 
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the topological trigger is a quick change of the global topology, which dictates the fast 
reconnection of collisional or collisionless origin. Thus the topological trigger is the most 
appropriate nomenclature to emphasize the basic nature of the topological effect involved, 
and it is a welcome usage. 

3.2 How does the topological trigger work? 

Let us trace how a rapid rearrangement of the field topology occurs under conditions of slow 
evolution of photospheric sources with the total charge equals zero. We shall arbitrary fix 
the positions of three charges, while we move the fourth one along an arbitrary trajectory 
in the plane Q. 

Let an initial position of the moving charge corresponds to the values of the topological 
indices Itop = — 1 and + 1 for the zeroth points X\ (Type A— with X z > 0) and X2 
(Type A+ with < 0) respectively (Figure [2]). Thus two points in the plane Q have 
different indices. Equation (115p is satisfied. The separator is the field line connecting 
these points without a coronal null. This field line emerges from the point X\ and is 
directed along the separator to the point X%. The value and sign of the magnetic-potential 
difference between the zeroth points are determined from the relation 



where the integral is taken along the separator. 1^2 — ^1 < for the field shown in Figure [2j 
According to Gorbachev et al. [33], the moving charge can arrive in a narrow region 
(let us call it the region TT) such that both points in the plane Q will have the same 
indices (Figure [TT]) . It follows from the general 3D Equation (|15|) that in this case there 
must also exist two zeroth points outside the plane. They are arranged symmetrically 
relative to the plane Q (the plane z = in this Section) . Figure [12] illustrates how these 
additional points appear. 

Before the start of trigger, the moving charge is outside of the region TT, the in- 
dex I top = — 1 and the eigenvalue A^ > at the non-degenerate zeroth point X\ (Fig- 
ure [12^,). When the moving charge crosses the boundary of the region TT, the eigen- 
value A z at the point X\ vanishes (Figure [12b): 



The point becomes degenerate. At this instant, another pair of zeroth points is born from 
the point X\ (Figure [12b) . We shall consider only one of them, the point X c in the upper 
half-space z > 0. This non-degenerate point travels along the separator and merges with 
the point X2 in the plane z = when the moving charge emerges from the region TT. At 
this instant, the eigenvalue \ z vanishes: 



2 




(29) 



1 



X z (X 1 ) = 0. 



(30) 



A 2 (A 2 ) = 0. 



(31) 
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Figure 12: Changes of the field pattern at the zeroth point X\. (a) An initial state is the 
non-degenerate point of the Type A— with A 2 > 0. (b) A degenerate hyperbolic point 
(line) with X z = at the beginning of trigger, (c) After the beginning of trigger, the 
pattern of field is the non-degenerate point of the Type A+ with \ z < and two zeroth 
points outside the plane z = 0. 

As a result of the process described, the direction of the field at the separator is reversed 
with the point A 2 of the Type A— with \ z > 0. After that, the moving charge is 
located outside the region TT, there are no zeroth points outside the plane z = 0. Thus 
Equations (f30j) and (f3Tj) determine the boundaries of the topological trigger region TT. 



y 




Figure 13: The topological trigger (nar- 
row dashed) region is restricted by two 
boundaries: curves C\ and C2 corre- 
sponding to Equations (f3"0"j) and (|3"Tj) . 



Figure 13 presents an illustrative sample of charge configuration. The fixed charges 
ei = — e2 = — e3 = 1 are located in the plane Q at the points (1; 0), (0; 0) and (0; 1), 
respectively. The charge e^ = 1 moves along a trajectory shown by a thin curve Tr. In the 
beginning of trigger, it crosses the boundary curve C\. Between the boundaries C\ and C2 
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there exists a curve C c at which ip2 — ipi = 0. The field at the separator changes sign at 
the point X c . Hence the sections X\X C and X C X2 of the separator make contributions 
of opposite signs to the integral (|29p . These contributions exactly compensate each other 
for a position of the point X c when the trajectory Tr crosses the curve C c . 

Typically the region TT is narrow. That is why small shifts of the moving charge 
within this region lead to large shifts of the zeroth point X c along the separator above the 
plane z = just creating a global bifurcation. Using the analogy with hydrodynamics, we 
see that the separatrix plane (y, z) in Figure fT2a. which plays the role of a "hard wall" for 
"flowing in" magnetic flux, is quickly replaced by the orthogonal "hard wall" , the separat- 
rix plane (x,z) in Figure [T2b. Thus the topological trigger drastically changes directions 
of magnetic fluxes in an AR as illustrated by Figure 14. 

Another specific sample is a charge arrangement along a straight line, e.g., the axis x 
in Figure [TJ In this case, owing to the axial symmetry, the entire separator consists of 
zeroth points. Thus we have a zeroth line. By using an inversion transformation (e.g. 
Landau et al. [38]): 




(32) 



where R is the inversion radius, it is easy to show that the axial symmetry is not a 
necessary condition for the appearance of zeroth line. Moreover, the new zeroth line also 
represents a circle centered in the plane z = 0. Thus the 3D zeroth lines of the magnetic 
field can exist if the sunspots do not lie in a straight line. 

Up to now, we have considered the travel of one charge while the coordinates and 
magnitudes of the other three charges were fixed. It is obvious, however, that all the 
foregoing remains in force in the more general case of variation of the charge configuration. 
It follows from the results that a slow evolution of the configuration of field sources in the 
photosphere can lead to a rapid rearrangement of the global topology in ARs in the corona. 
The phenomenon of topological trigger is necessary to model the large eruptive flares. 

4 Basic physics of reconnection in the corona 
4.1 Super-hot turbulent-current layers 

Coulomb collisions do not play any role in the super-hot turbulent- current layers (SHTCL). 
So the plasma inside the SHTCL has to be considered as collisionless (see Somov [19]). 
The concept of an anomalous resistivity, which originates from wave-particle interactions, 
is then useful to describe the fast conversion from field energy to particle energy. Some of 
the general properties of such a collisionless reconnection can be examined in a frame 
of a self-consistent model which makes it possible to estimate the main parameters of the 
SHTCL. 

As it reconnects, every magnetic-field line at first penetrates into the current layer with 
velocity vq and later on it moves out with velocity vi as illustrated by Figure [15j Basing 
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(a) 




Figure 14: Global changes 
of the magnetic field con- 
figuration related to the 
topological trigger effect in 
the model of an active re- 
gion as illustrated by Fig- 
ure 13. (a) An initial state 
before the trigger; the ze- 
roth point X\ is the non- 
degenerate point with a ver- 
tical component of mag- 
netic field B z > 0. (b) Af- 
ter beginning of the trig- 
ger, the point X\ is the 
non-degenerate point with 
B z < 0; there are two 
zeroth points outside the 
plane z = 0, one of them 
is shown as X c above the 
plane z = 0. (c) The end of 
the topological trigger pro- 
cess. 
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Figure 15: A magnetically non-neutral reconnecting layer: the electric current distribution 
is schematically shown by the shadow, the dotted boundary indicates the field lines going 
through the current layer. 2a is the thickness of the current layer, 26 is its width. 2a out 
is an effective cross-section for the outflows of energy and mass along the reconnected 
magnetic-field lines. 



on the mass, momentum and energy conservation laws, we write the following relations: 

n v b = n s v 1 a out , (33) 

2n k B T + = n s k B T (l + ^ , (34) 



1\ 1 
2 



n s k B T 1 + - ) = -Mn s v? + 2n k B T , (35) 



Xei£Z g + £iZ,e=£%!e+Cr, (36) 

(1 - Xef ) SZg + £ t,i = £Z + K?* ■ ( 37 ) 

Here no and n s are the plasma densities outside and inside the layer. To is the temperature 
of inflowing plasma outside the layer, T = T e is an effective electron temperature inside 
it, the ratio 6 = T e /Ti, Tj is an effective temperature of ions. The velocity of electric drift 
of plasma in direction to the layer 

v = V d = c^, (38) 

-DO 

and 

vi = V A -= , B ° (39) 
is the velocity of the plasma outflow. 
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The continuity Equation (|33|) and the energy Equations ([36]) and ([37|) are of integral 
form for a quarter of the layer assumed to be symmetrical and for a unit length along the 
current. The left-hand sides of (|36p and (|37p contain the magnetic energy flux 

£ ™ g = ^vob, (40) 

which creates heating of electrons and ions due to their interactions with waves. A relative 
fraction Xei of the heating is consumed by electrons, while the remaining fraction (1 — Xef) 
goes to the ions. 

The electron and ion temperatures of the plasma inflowing to the layer are the same. 
Hence, the fluxes of the electron and ion thermal energies are also the same: 

£lZ, e = £^ = ln k B T -v b. (41) 

The factor 5/2 appears because, for the particles of kind k (electrons, protons, other ions), 
the mean kinetic energy flux of chaotic motion, n k m k v k w k , where the heat function per 
unit mass or, more exactly, the specific enthalpy is 

p k 5 k B T k 

w k = e k -\ = • 

Pk 2 m k 

Because of the difference between the temperatures of electrons and ions in the out- 
flowing plasma, the electron and ion thermal energy outflows differ: 



??;,», =\n s k B T. Vl a° u \ £Z = l n sk B ^- Vl a^ . (12) 



The ion kinetic energy flux from the layer 

K,° ut = i Mn s v? ■ Vl a out (43) 

is important in the energy balance (j37|) . As to the electron kinetic energy, it is negligible 
and disregarded in (|36p . However, electrons play the dominant role in the heat conductive 
cooling of a SHTCL: 

cf = f M (0) ns{k ^ 2 a ° ut - ( 44 ) 

Here f M (8) is the function which allows us to consider the field-aligned anomalous thermal 
flux depending on the the ratio 9. 

Under the conditions derived from the Yohkoh and RHESSI data, contributions to the 
energy balance are not made either by the energy exchange between the electrons and the 
ions due to collisions, the thermal flux across the magnetic field, and the energy losses for 
radiation. The field-aligned thermal flux becomes anomalous and plays the dominant role 
in the cooling of electron component inside the SHTCL. All these properties are typical 
for collisionless l super-hoV (T e > 30 MK) plasma. 
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Under the same conditions, the anomalous conductivity cr e f in the Maxwell equation 
for curlB 

f^ = a ci E , (45) 

as well as the relative fraction Xef of the direct heating consumed by electrons, are deter- 
mined by the wave-particle interaction inside the SHTCL and depend on a type of plasma 
turbulence and its regime. 



4.2 Plasma turbulence inside the SHTCL 

In Equation (|4"5]) it is convenient to replace the effective conductivity cr e f by effective 
resistivity r] e f. 

^ = — • (46) 

47TO 7/ ef 

In general, the partial contributions to the resistivity may be made simultaneously by 
several processes of electron scattering by different sorts of waves, so that 

Vef = J2r)k- (47) 

k 

The relative share of the electron heating Xef is also presented as a sum of the respec- 
tive shares Xk of the feasible processes taken, of course, with the weight factors r?fc/7? e f 
which defines the relative contribution from one or another process to the total heating of 
electrons inside the SHTCL: 

Xcf = E — Xk ■ (48) 

T ^ 

In usual practice (e.g. Somov [26]), the sums (|47p and ([48 p consist of no more than 
two terms, either of which corresponds to one of the turbulent types or states. 



4.3 Marginal and saturation regimes 

When the electron current velocity u = j/n s exceeds a critical value, the instabilities due 
to current flow of electrons appear. A rapid decrease in the plasma conductivity occurs 
and anomalous resistivity arises (Kadomtsev [39], Artsimovich and Sagdeev [40]). The 
condition needed for current instability in a RCL is that the layer thickness is of the 
order of the ion gyroradius (Syrovatskii [17]). Thus the turbulent current layers must be 
sufficiently thin. According to Syrovatskii [41], the development of a thin current layers 
(TCL) in the solar atmosphere leads to plasma turbulence and correspondingly to a fast 
rate of field dissipation in flares with the heating of plasma to high temperatures, the 
high-velocity plasma ejections, and the acceleration of particles to high energies. 

Let us consider the turbulence which is due only to the ion-cyclotron (ic) and ion- 
acoustic (ia) instabilities. The sums (j47p and (|48|) consist of no more than two terms; either 
of which corresponds to one of the said turbulence types. In reality, the ion-cyclotron waves 
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Figure 16: The relative velocity of electrons, U = u/Vr e as a function of the ratio 8 = T e /Ti 
inside the SHTCL. The individual 'arcs' of the curves correspond to four different regimes 
of turbulence. The dashed continuations of the 'arcs' are the critical velocities of excitation 
of corresponding waves. The arrow A shows that, in equilibrium plasma with 9 = 1, the 
ion-cyclotron waves (the curved section KO) are excited earlier than the ion-acoustic waves 
(the arc LO). 

In the case of the marginal regime of turbulence, the velocity of electrons 

u = (49) 
en s r] ei 

coincides with the critical velocity Uk of the A:-type wave excitation. Hence 

U k (6) = U k mar (*) = ^ ■ (50) 

For example, the ion-cyclotron instability becomes enhanced when the velocity u is not 
lower than the critical value Uj c of the ion-cyclotron waves. In the marginal regime 

ujr (o) = ^. (51) 

This function is shown as the curved section KO in Figure [TBI 

As long as the ion-cyclotron waves are not saturated, the velocity u remains approx- 
imately equal to Uj c and thus it is possible to calculate the resistivity 77 e f from Equa- 
tion (149p . Therefore, if the velocity u in the SHTCL does not exceed the ion-acoustic 
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wave excitation threshold U{ a , only the ion-cyclotron waves will contribute to anomalous 
resistivity 77 e f and the factor Xet- 

In the marginal regime, the wave-particle interaction is quasilinear. However, in the 
case of sufficiently strong electric field Eq, the nonlinear interactions become important, 
thereby giving rise to another state of turbulence. In this regime, Ohm's law can no 
longer define the resistivity, but determines the the electron current velocity. Regarding 
the resistivity, this is inferred from the turbulence saturation level. In the saturated regime, 
Uf, (9) must be replaced by certain functions U^ Lt (9), shown by the curved section KL 
in Figure [TUl and U i s ( ^ lt (9), shown by the arcs MQ for the ion-cyclotron and ion-acoustic 
turbulence, respectively. 

Calculations show that the energy release power and the reconnection rate, which are 
necessary for solar flares to be accounted for, can be obtained in the marginal regime of 
ion-acoustic turbulence (see Somov [26,19]). 

5 The collapsing magnetic trap effects 
5.1 Fast and slow reconnection 

Collapsing magnetic traps are formed by the process of collisionless reconnection in the 
solar atmosphere (Somov and Kosugi [42]). Figure [TTJ illustrates two possibilities. Fast 
(Figure fTTh ) and slow (Figure fl~7b) modes of reconnection are sketchy shown in the corona 
above the magnetic obstacle, the region of a strong magnetic field, which is observed in 
SXRs as a flare loop (shaded). 

In the first case, let us assume that both feet of a reconnected field loop path through 
the shock front (SW in Figure fTTh ) ahead the obstacle. Depending on the velocity and 
pitch-angle, some of the particles preaccelerated by the SHTCL may penetrate through the 
magnetic-field jump related to the shock or may be reflected. For the particles reflected 
by the shock, the magnetic loop represents a trap whose length L(t), the distance between 
two mirroring points at the shock front, measured along a magnetic-field line, decreases 
from its initial value L(0) ~ 2Lq to zero (the top of the loop goes through the shock front) 
with the velocity ~ 2v\. Therefore, the lifetime of each collapsing trap t± « Lq/v\. 

In the case of slow reconnection, there is no a shock wave, and the trap length L{t) 
is the distance between two mirroring points (Mi and M2 in Figure [TTb). measured along 
a reconnected magnetic-field line. In both cases, the electrons and ions are captured in 
a trap whose length decreases. So the particles gain energy from the increase in parallel 
momentum. 

Thus, in the first approximation, we neglect collisions of particles ahead of the shock 
wave (Figure [TTa) or in the trap without a shock (Figure [TTb). In both cases, the particle 
acceleration can be demonstrated in a simple model - a long trap with short mirrors 
(Figure [T8|) . The decreasing length L(t) of the trap is much larger than the length l m of 
the mirrors; the magnetic field B = Bi is uniform inside the trap but grows from Bi to 
B2 in the mirrors. The quantity B2/B1 is called the mirror ratio; the larger this ratio, 
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Figure 17: Plasma flows related to a super-hot turbulent-current layer (SHTCL): the 
inflows with a relatively low velocity vo, the downward outflow with a super- Alfven ve- 
locity vi. (a) SW is the shock wave above the magnetic obstacle. V2 is the post-shock 
velocity, vy is the velocity of spreading of the compressed plasma along the field lines 
toward the feet of the loop, (b) The supra-arcade downflow and collapsing trap without 
a shock. Mi and M2 are the mirroring points where the field becomes sufficiently strong 
to reflect fast particles above the chromosphere (Ch). 
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the higher the particle confinement in the trap. The validity conditions for the model are 
discussed by Somov and Bogachev [43]. 



5.2 The first-order Fermi-type acceleration 

We consider the traps for those the length scale and timescale are both much larger than 
the gyroradius and gyroperiod of an accelerated particle. Due to strong separation of 
length and timescales, the magnetic field inside the trap can be considered as uniform 
and constant (for more detail see Somov and Bogachev [43]). If so, then the longitu- 
dinal momentum of a particle increases with a decreasing length L(t), in the adiabatic 
approximation, as 

vS) = P ~f- (52) 

Here I = L(t)/L(0) is the dimensionless length of the trap. The transverse momentum is 
constant inside the trap, 

P±=P±o, (53) 
because the first adiabatic invariant is conserved: 

v 2 

^ = const . (54) 
Thus the kinetic energy of the particle increases as 

The time of particle escape from the trap, I = l es , depends on the initial pitch- angle 6$ 
of the particle and is determined by the condition 

tg0 O = — <^-, (56) 

P\\o nl es 




Figure 18: Two main effects in a collapsing trap, (a) Magnetic mirrors move toward each 
other with velocity v m . (b) Compression of the trap with velocity v t . 
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where 




*-(£-f ■ 

at the time of its escape 



(57) 



2m V y 2m Bi 



(58) 



One can try to obtain the same canonical result by using more complicated approaches. For 
example, Giuliani et al. [44] numerically solved the drift equations of motion. However it is 
worthwhile to explore first the simple analytical approach to investigate the particle energi- 
zation processes in collapsing traps in more detail before starting to use more sophisticated 
methods and simulations. 

5.3 The betatron acceleration in a collapsing trap 

If the thickness of the trap decreases with its decreasing length, then the strength of the 
field Bi inside the trap increases as a function of I, say B\(l). In this case, according 
to (|54p . the transverse momentum increases simultaneously with the longitudinal momen- 
tum (1521) : 



Here B\ = B\(l) is the initial (at / = 1) value of the field inside the trap. The kinetic 
energy of a particle 



increases faster than that in the absence of trap contraction, see (|55p . Therefore it is 
natural to assume that the acceleration efficiency in a collapsing trap also increases. 
However, as the trap is compressed, the loss cone becomes larger (Figure \19\i . 



Consequently, the particle escapes from the trap earlier. 

On the other hand, the momentum of the particle at the time of its escape satisfies 
the condition 




(59) 




(60) 




(61) 



Pll (l) = R(l) P± (l), 



(62) 



where 




(63) 
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Figure 19: The betatron effect in a collapsing magnetic trap. As the trap is compressed 
with velocity Vj, the loss cone becomes larger. A particle escapes from the trap earlier 
with an additional energy due to betatron acceleration. 



Hence, using (159D . we determine the energy of the particle at the time of its escape from 
the trap 

P^l) 2 / Dm2 , A Plo Bi(l) B 2 p 2 n B, 



K P 



2 m 



(my 



+1 



2m Si Bi(Z) ~ 2m Si ' 1 J 

The kinetic energy (I64p . that the particle gains in a collapsing trap with compression, 
is equal to the energy (|58p in a collapsing trap without compression, i.e. without the 
betatron effect. 

Thus the compression of a collapsing trap (as well as its expansion or the transverse 
oscillations) does not affect the final energy that the particle acquires during its accelera- 
tion. The faster gain in energy is exactly offset by the earlier escape of the particle from 
the trap (Somov and Bogachev [43]). The acceleration efficiency, which is defined as the 
ratio of the final (I = li s ) and initial (/ = 1) energies, i.e. 



K es 

K(l) 



Plo 



P 2 ±o+P 2 o 



B2 

B! 



P±o 



Eh 
Bi 



(65) 



depends only on the initial mirror ratio B2/B1 and the initial particle momentum or, to 
be more precise, on the ratio p j_ /po- The acceleration efficiency (f65j) does not depend on 
the compression of collapsing trap and the pattern of decrease in the trap length either. 

It is important that the acceleration time in a collapsing magnetic trap with compres- 
sion can be much shorter than that in a collapsing trap without compression. For example, 
if the cross-section area S(l) of the trap decreases proportionally to its length /: 



S(l) = S(l)l, 
then the magnetic field inside the trap 

B 1 (l) = B 1 (l)/l, 



(66) 



(67) 
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and the effective parameter 



1 - A 1 / 2 

i?(Z) = ( i? 2 — ) , (68) 



where R is define by formula (|57p . At the critical length 

the magnetic field inside the trap becomes equal the field in the mirrors, and the magnetic 
reflection ceases to work. If, for certainty, B2/B1 = 4, then l cr = 1/4. So contraction of 
the collapsing trap does not change the energy of the escaping particles but this energy is 
reached at an earlier stage of the magnetic collapse when the trap length is finite. In this 
sense, the betatron effect increases the actual efficiency of the main process - the particle 
acceleration on the converging magnetic mirrors. 



5.4 The betatron acceleration in a shockless trap 

If we ignore the betatron effect in a shockless collapsing trap, show in Figure \T7b . then 
the longitudinal momentum of a particle is defined by the formula (instead of (|52p ) 

p I, (t) W p „ (0) A* + h) p || (0) ffl±M , when t^t x . (70) 

The particle acceleration on the magnetic mirrors stops at the time t\ = l\jv\ at a finite 
longitudinal momentum that corresponds to a residual length (J2 in Figure 117b) of the 
trap. 

Given the betatron acceleration due to compression of the trap, the particle acquires 
the same energy (|58p by this time or earlier if the residual length of the trap is comparable 
to a critical length l cr determined by a compression law (Somov and Bogachev [43]). 
Thus the acceleration in shockless collapsing traps with a residual length becomes more 
plausible. The possible observational manifestations of such traps in the X-ray and optical 
radiation are discussed by Somov and Bogachev [43]. The most sensitive tool to study 
behavior of the electron acceleration in the collapsing trap is radio radiation. We assume 
that wave-particle interactions are important and that two kinds of interactions should be 
considered in the collapsing trap model. 

The first one is resonant scattering of the trapped electrons, including the loss-cone 
instabilities and related kinetic processes (e.g. Benz [45], Chapter 8). Resonant scattering 
is most likely to enhance the rate of precipitation of the electrons with energy higher that 
hundred keV, generating microwave bursts. The lose-cone instabilities of trapped mildly- 
relativistic electrons (with account taken of the fact that there exist many collapsing field 
lines at the same time, each line with its proper time-dependent loss cone) would provide 
excitation of waves with a very wide continuum spectrum. In a flare with a slowly-moving 
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upward coronal HXR source, an ensemble of the collapsing field lines with accelerated 
electrons would presumably be observed as a slowly moving type IV burst with a very 
high brightness temperatures and with a possibly significant time delay relative to the 
chromospheric footpoint emission. 

The second kind of wave-particle interactions in the collapsing trap-plus-precipitation 
model is the streaming instabilities (including the current instabilities related to a return 
current) associated with the precipitating electrons. 

5.5 Some observational results 

In order to interpret the temporal and spectral evolution and spatial distribution of HXRs 
in flares, a two-step acceleration was proposed by Somov and Kosugi [42] with the second 
step of acceleration via the collapsing magnetic-field lines. The Yohkoh HXT observations 
of the Bastille-day flare (Masuda et al. [46]) clearly show that, with increasing energy, 
the HXR emitting region gradually changes from a large diffuse source, which is located 
presumably above the ridge of soft X-ray arcade, to a two-ribbon structure at the loop 
footpoints. This result suggests that electrons are in fact accelerated in the large system 
of the coronal loops, not merely in a particular one. This seems to be consistent with the 
RHESSI observations of large coronal HXR sources; see, for example, the X4.8 flare of 
2002 July 23 (see Figure 3 in Lin et al. [47]). 

Efficient trapping and continuous acceleration also produce the large flux and time 
lags of microwaves that are likely emitted by electrons with higher energies, several hun- 
dred keV (Kosugi et al. [48]). The lose-cone instabilities (Benz [45]) of trapped mildly- 
relativistic electrons in the system of many collapsing field lines (each line with its proper 
time-dependent lose cone) presumably can provide excitation of radio-wave with a very 
wide continuum spectrum. 

Qiu et al. [49] presented a comprehensive study of the X5.6 flare on 2001 April 6. 
Evolution of HXRs and microwaves during the gradual phase in this flare exhibits a sep- 
aration motion between two footpoints, which reflects the progressive reconnection. The 
gradual HXRs have a harder and hardening spectrum compared with the impulsive com- 
ponent. The gradual component is also a microwave-rich event lagging the HXRs by tens 
of seconds. The authors propose that the collapsing-trap effect is a viable mechanism that 
continuously accelerates electrons in a low-density trap before they precipitate into the 
footpoints. 

Imaging radio observations (e.g. Li and Gan [50]) should provide another way to 
investigate properties of collapsing magnetic traps. It is not simple, however, to understand 
the observed phenomena relative to the results foreseen by theory. With the incessant 
progress of magnetic reconnection, the loop system newly formed after reconnection will 
grow up, while every specific loop will shrink. Just because of such a global growth of 
flare loops, it is rather difficult to observe the downward motion of newly formed loops. 
The observations of radio loops by Nobeyama Radiohcliograph (NoRH) are not sufficient 
to resolve specific loops. What is observed is the whole region, i.e., the entire loop or 
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the loop top above it. Anyway, combined microwave and HXR imaging observations are 
essential in the future. 

6 Open issues of reconnection in flares 

The existing models of reconnection in the solar corona can be classified in two wide 
groups: global and local ones (Figure 1201) . The global models are used to describe actual 
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Figure 20: Models of magnetic reconnection in the solar atmosphere. 

ARs or even complexes of activity in different approximations and with different accuracies 
(Somov [51], Gorbachev and Somov [7], Bagala et al. [52], Antiochos [31], Aschwanden et 
al. [53], Morita et al. [54], Longcope et al. [55], Barnes [34], Longcope and Beveridge [56], 
Ugarte-Urra et al. [35]). We make no attempt to review all the models but just remark 
that the main advantage of the global models is direct comparison between the results of 
computation and the observed large-scale patterns. For example, the 1 rainbow reconnec- 
tion 1 model (Somov [57]) is used to reproduce the main features of magnetic fields in the 
corona, related to the large-scale vortex flows in the photosphere. 

The advantage of the local models is that they take kinetic effects into account and 
allow us to develop the basic physics of reconnection in flares. In general, many analytical, 
numerical, and combined models of reconnection exist in different approximations and with 
different levels of self-consistency (see Biskamp [15], Priest and Forbes [58]). It becomes 
more and more obvious that collisionless reconnection in a rarefied plasma is a key process 
in flares. This process was introduced by Syrovatskii [59] as a dynamic dissipation of field 
in a current layer and leads to fast conversion from field energy to particle energy, as well 
as a topological change of the magnetic field (e.g. Horiuchi et al. [60], Yamada et al. [16]). 

General properties and parameters of the collisionless reconnection can be examined 
in a frame of models based on the mass, momentum, and energy conservation laws (Sec- 
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tion 2]) . A particular feature of the models is that electrons and ions are heated by wave- 
particle interactions in a different way. The magnetic-field-aligned thermal flux becomes 
anomalous and plays the role in the cooling of electrons in a super-hot turbulent- current 
layer (SHTCL). These properties seem to be typical for the conditions derived from the 
observations by Yohkoh and RHESSI (e.g. Joshi et al. [61], Liu et al. [62]), the SOXS 
mission (Jain et al. [63]). Unfortunately, the local models are not incorporated in the 
global consideration of reconnection in the corona. Only a few first steps have been made 
in this direction (see Somov [19]). 

Spacecraft observations of collisionless reconnection in the magnetotail and magne- 
topause as well as recent simulations have shown the existence of thin current layers 
(TCLs), with scale lengths of the order a few electron skin depth. In the electron MHD 
model of such TCLs, reconnection is facilitated by electron inertia which breaks the frozen- 
in condition (Jain and Sharma [64]). The simulations demonstrate the instability of the 
whistler-like mode which presumably plays the role of the ion-cyclotron instability in the 
SHTCL model. Stark et al. [65] combine the TCL models with a global model of magnetic 
field in the magnetosphere. 

As for particle acceleration in flares, it goes in two steps. During the first one, the 
ions and electrons are accelerated by the DC electric field in a super-hot turbulent-current 
layer (Litvinenko and Somov [66,67]). During the second step, the energy of particles 
trapped in collapsing traps additionally increases by the first-order Fermi mechanism (So- 
mov and Kosugi [42]) and the betatron mechanism (Somov and Bogachev [43], Karlicky 
and Kosugi [68]). In a collisionless approximation, Bogachev and Somov [69] show that 
the collapsing trap model predicts two types of coronal HXR sources: thermal and non- 
thermal. Thermal sources are formed in traps dominated by the betatron mechanism. 
Non-thermal sources with power-law spectra appear when electrons are accelerated by the 
Fermi mechanism. These results are of interest in interpreting the RHESSI observations. 

Future models of solar flares should join global and local properties of reconnection 
under coronal conditions. For example, chains of plasma instabilities, including kinetic 
instabilities, can be important for our understanding the types and regimes of plasma 
turbulence inside the collisionless current layers with a longitudinal magnetic field. In 
particular it is necessary to evaluate better the anomalous resistivity and selective heat- 
ing of particles in such a SHTCL. Heat conduction is also anomalous in the super-hot 
plasma. Self-consistent solutions of the reconnection problem will allow us to explain the 
energy release in flares, including the open question of the mechanism or combination of 
mechanisms which explains the observed acceleration of electrons and ions to high energy. 

To understand the 3D structure of reconnection in flares is one of the most urgent 
problems. Actual flares are 3D dynamic phenomenon of electromagnetic origin in a highly- 
conducting plasma with a strong magnetic field. RHESSI and Hinode observations offer 
us the means to check whether phenomena predicted by topological models (such as the 
topological trigger) do occur [70]. However some puzzling discrepancies may also exist, 
and further development of realistic 3D models is required. 
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